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(Dated: September 25, 2006) 

We present a measurement of the fraction /+ of right-handed W bosons produced in top quark decays, based 
on a candidate sample of ti events in the i?+jets and dilepton decay channels corresponding to an integrated 
luminosity of 370 pb ~ collected by the D0 detector at the Fermilab Tevatron pp Collider at y/s = 1.96 Te V. We 
reconstruct the decay angle 9* for each lepton. By comparing the cos 6* distribution from the data with those for 
the expected background and signal for various values of /+ , we find /+ = 0.056±0.080 (stat) ±0.057 (syst). 
(/+ < 0.23 at 95% C.L.), consistent with the standard model prediction of /+ = 3.6 x 10~ 4 . 



PACS numbers: 14.65 .Ha, 14.70.Fm, 12.15.Ji, 12.38.Qk, 13.38.Be, 13.88.+e 
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The top quark is by far the heaviest of the known fermions 
and is the only one that has a Yukawa coupling of order unity 
to the Higgs boson in the standard model. We search for ev- 
idence of new physics in t — > Wb decay by measuring the 
helicity of the W boson. In the standard model, the top quark 
decays via the V — A charged current interaction, almost al- 
ways to a W boson and a b quark. For any linear combination 
of V and A currents at the t — > Wb vertex, the fraction /o of 
longitudinally-polarized W bosons is 0.697 ± 0.012 JJJ] at the 
world average top quark mass m t of 172.5 ± 2.3 GeV |2]. 

In this analysis, we fix /o at 0.70 and measure the positive 
helicity fraction / + . In the standard model, / + is predicted at 
next-to-leading order to be 3.6 x 10 -4 yj. A measurement 
of /_)_ that differs significantly from this value would be an 
unambiguous indication of new physics. For example, an / + 
value of 0.30 would indicate a purely V + A charged current 
interaction. 

Measurements of the b — ► decay rate have indirectly 
limited the V + A contribution in top quark decays to less 
than a few percent |4]. Direct measurements of the V + A 
contribution are still necessary because the limit from b — > sj 
assumes that the electroweak penguin contribution is domi- 
nant. Direct measurements of the longitudinal fraction found 
fa = 0.91 ± 0.39 |5] and f = 0.56 ± 0.31 H. Di- 
rect measurements of /+ have set limits of / + < 0.18 0], 
/+ < 0.24 H, and /+ < 0.25 |9] at the 95% C.L. The an- 
alysis presented in this Letter improves upon that reported in 
Ref. 1 9] by using a larger data set, including the dilepton de- 
cay channel of the ti pair, and employing enhanced analysis 
techniques. 

The angular distribution of the down-type decay products of 
the W boson (charged lepton or d, s quark) in the rest frame 
of the W boson can be described by introducing the decay an- 
gle 8* of the down-type particle with respect to the top quark 
direction. The dependence of the distribution of cos 8* on / + , 

c J (c e ,)«2(l- C 2,)/o + (l-c e 2 /- +(l + c e ,) 2 /+ (1) 

where c* g — cos 8*, forms the basis for our measurement. We 
proceed by selecting a data sample enriched in ti events, re- 
constructing the four vectors of the two top quarks and their 
decay products, and then calculating cos 8* . This distribution 
in cos 8* is compared with templates for different / + values, 
suitably corrected for background and reconstruction effects, 
using a binned maximum likelihood method. In the ^+jets 
channel, the kinematic reconstruction is done with a fit that 
constrains the W boson mass to its measured value and the 
top quark mass to 175 GeV, while in the dilepton channel, the 
kinematics are solved algebraically with the top quark mass 
fixed to 172.5 GeV. 

The D0 detector fioll comprises three main systems: the 
central tracking system, the calorimeters, and the muon sys- 
tem. The central-tracking system is located within a 2 T 
solenoidal magnet. The next layer of detection involves three 
liquid-argon/uranium calorimeters: a central section covering 
pseudorapidities \rj\ < 1, and two end calorimeters ex- 
tending coverage to \rj\ w 4, all housed in separate cryostats. 



The muon system is located outside the calorimetry, and con- 
sists of a layer of tracking detectors and scintillation trigger 
counters before 1.8 T toroids, followed by two similar layers 
after the toroids. 

This measurement uses a data sample recorded with the 
D0 experiment and corresponds to an integrated luminosity 
of about 370 pb _1 of pp collisions at y/s = 1.96 TeV. The 
data sample consists of ti candidate events from the i?+jets 
decay channel ti — > W + W~bb — > tvqq'bb and the dilepton 
channel ti — > W + W~bb — » tvl'v'bb, where I and I' are elec- 
trons or muons. The ^+jets final state is characterized by one 
charged lepton, at least four jets (two of which are b jets), and 
significant missing transverse energy ($t)- The dilepton final 
state is characterized by two charged leptons of opposite sign, 
at least two jets, and significant fSr- 

We simulate ti signal events with mt = 172.5 GeV for 
different values of /+ with the ALPGEN Monte Carlo (MC) 
program 112)1 for the parton-level process (leading order) and 
PYTHIA 11311 for gluon radiation and subsequent hadroniza- 
tion. As the interference term between V — A and V + A is 
suppressed by the small mass of the b quark and is therefore 
negligible 1 14], samples with f + = 0.00 and f + = 0.30 are 
used to create cos 8* templates for any / + value by a linear 
interpolation of the templates. 

The MC samples used to model background events with 
real leptons are also generated using ALPGEN and PYTHIA. 
Backgrounds in the £+jets channel arise predominantly from 
W^+jets production and multijet production where one of the 
jets is misidentified as a lepton and spurious $t appears due 
to mismeasurement of the transverse energy in the event. 

The £+jets event selection fl5ll requires an isolated lepton 
(e or /i) with transverse momentum px > 20 GeV, no other 
lepton with p T > 15 GeV in the event, $ T > 20 GeV, and 
at least four jets. Electrons are required to have < 1.1 and 
are identified by their energy deposition and isolation in the 
calorimeter, their transverse and longitudinal shower shapes, 
and information from the tracking system. Also, a discrimi- 
nant combining the above information must be consistent with 
the expectation for a high-pT isolated electron 1 15]. Muons 
are identified using information from the muon and tracking 
systems, and must satisfy isolation requirements based on the 
energies of calorimeter clusters and the momenta of tracks 
around the muon. They are required to have |?7| < 2.0 and 
to be isolated from jets. Jets are reconstructed using the Run 
II mid-point cone algorithm with cone radius 0.5 I1 1611 . and are 
required to have rapidity \y\ < 2.5 and p? > 20 GeV. 

We determine the number of multijet background events 
N m j from the data, using the technique described in Ref. fl5ll . 
We calculate N m j for each bin in the cos 8* distribution from 
the data sample to obtain the multijet cos 8* templates. 

To discriminate between ti pair production and back- 
ground, a discriminant T> with values in the range to 1 is 
calculated using input variables which exploit differences in 
kinematics and jet flavor. The kinematic variables considered 
are: Ht (defined as the scalar sum of the jet pt values), the 
minimum dijet mass of the jet pairs rrijj m i n , the \ 2 from the 
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kinematic fit, the difference in azimuthal angle Acj) between 
the lepton and $t directions, and aplanarity A and sphericity 
S 11711 (calculated from the four leading jets and the lepton). 
Only the four leading jets in p-p are considered in computing 
these variables. 

We utilize the fact that background jets arise mostly from 
light quarks or gluons while two of the jets in tt events arise 
from b quarks by considering the impact parameters with re- 
spect to the primary vertex of all tracks within the jet cone. 
Based on these values, we calculate the probability Ppy for 
each jet to originate from the primary vertex. We then average 
the two smallest Ppy values to form a continuous variable 
(Ppv) that tends to be small for tt events and large for back- 
grounds. This approach results in similar background dis- 
crimination but better efficiency than applying a simple cut 
on Ppy. 

The discriminant is built separately for the e+jets and 
/i+jets channels, using the method described in Refs. | la Jlal . 
Background events tend to have T> values near 0, while tt 
events tend to have values near 1. We consider all possible 
combinations of the above variables for use in the discrimi- 
nant, and all possible requirements on the T> value, and choose 
the variables and T> criterion that give the smallest expected 
uncertainty on /+. In the e+jets channel, S, Hp, (Ppy), and 
X 2 are used, and T> is required to be > 0.65. In the /i+jets 
channel, A, Hp, nijjmin, (Ppy), \ 2 , an d Acf> are used, and T> 
is required to be > 0.80. In both channels the efficiency for 
tt events to satistfy the T> requirement is independent of the 
value of /_)_ . 

We then perform a binned Poisson maximum likelihood fit 
to compare the observed distribution of events in T> to the sum 
of the distributions expected from tt, VK+jets, and multijet 
events. N m j is constrained to the expected value within the 
known uncertainty. The likelihood is then maximized with re- 
spect to the numbers of tt, W+jets, and multijet events, which 
are multiplied by the appropriate efficiency for the T> selection 
to determine the composition of the sample used for measur- 
ing cos 9*. 

In the dilepton channel, backgrounds arise from processes 
such as VFW^+jets or Z+jets. These processes are either rare 
or require false Jfip from mismeasurement of jet and lepton 
energy, allowing a good signal to background ratio to be at- 
tained using only kinematic selection criteria. The selection 
is detailed in Ref . 1 1 911 . Events are required to have two lep- 
tons with opposite charge and pp > 15 GeV and two or more 
jets with pt > 20 GeV and \y\ < 2.5. Additional criteria are 
applied in the ee and fj.fi channels to suppress Z — > It, and 
in the e/z channel the sum of the two leading jet pp's and the 
leading lepton pp must be greater than 122 GeV. We place a 
more stringent requirement on electron identification than is 
used in Ref. d. 

Table|l]lists the composition of each sample as well as the 
number of observed events in the data. We observe a dis- 
parity between the number of tt events in the e+jets chan- 
nel and /i+jets channel, which is unexpected since the selec- 
tion efficiencies for the two channels are similar. The statis- 



TABLE I: Number of events observed in each ti decay channel, the 
background level as determined by a fit to the T> distribution in the 
^+jets channels and the expectation from the background produc- 
tion rate and selection efficiency in the dilepton channels, and the 
expected signal yield assuming standard model tt production with a 
top quark mass of 175 GeV. 



Observed Background Expected tt 


e + jets 


51 


5.3 ±0.9 


32.9 


\i + jets 


19 


3.3 + 0.4 


26.4 


e/i 


15 


2.2 + 0.6 


8.9 


ee 


4 


0.8 + 0.2 


3.3 


tip, 


1 


0.4 + 0.1 


2.4 



tical significance of the discrepancy in the event distribution 
is slightly above 2a. The disparity appears to be a feature of 
the data sample used in this analysis, as it occurs regardless 
of the choice of variables used to define T>. Further, it has 
no direct impact on this analysis, which relies only upon the 
distribution of events in cosO* . 

The top quark and W boson four-momenta in the selected 
£+jets events are reconstructed using a kinematic fit which is 
subject to these constraints: two jets must form the invariant 
mass of the W boson, the lepton and the ftp together with the 
neutrino p z component must form the invariant mass of the 
W boson, and the masses of the two reconstructed top quarks 
must be 175 GeV. Among the twelve possible jet combina- 
tions, the solution with the minimal \ 2 from the kinematic fit 
is chosen; MC studies show this yields the correct solution in 
about 60% of all cases. The cos 9* distribution obtained in the 
£+jets data after the full selection and compared to standard 
and V + A model expectations is shown in Fig. [Ha). 

Dilepton events are rarer than £+jets events, but have the 
advantage that cos 9* can be calculated for each lepton, thus 
providing two measurements per event. The presence of two 
neutrinos in the dilepton final state makes the system kine- 
matically underconstrained. However, if a top quark mass is 
assumed, the kinematics can be solved algebraically with a 
four-fold ambiguity in addition to the two-fold ambiguity in 
pairing jets with leptons. For each lepton, we calculate the 
value of cos 9* resulting from each solution with each of the 
two leading jets associated with the lepton. To account for 
detector resolution we repeat the above procedure 100 times, 
fluctuating the jet and lepton energies within their resolutions 
for each iteration. The average of these values is taken as 
the cos 9* for that lepton. The cos 9* distribution obtained in 
dilepton data is shown in Fig.Qb). 

We compute the binned Poisson likelihood L(f + ) for the 
data to be consistent with the sum of signal and background 
templates at each of seven chosen / + values. A parabola is 
fit to the — In [L (/+)] points to determine the likelihood as a 
function of / + . 

Systematic uncertainties are evaluated in ensemble tests by 
varying the parameters (see Table |ll|l which can affect the 
shapes of the cos 9* distributions or the relative contribution 
from signal and background sources. Ensembles are formed 
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FIG. 1: cos 8* distribution observed in (a) ^+jets and (b) dilepton 
events. The standard model prediction is shown as the solid line, 
while a model with a pure V + A interaction would result in the 
distribution given by the dashed line. 



by drawing events from a model with the parameter under 
study varied. These are compared to the standard cos 9* tem- 
plates in a maximum likelihood fit. The average shift in the 
resulting / + value is taken as the systematic uncertainty and 
is shown in Table ITTI The total systematic uncertainty is then 
taken into account in the likelihood by convoluting the lat- 
ter with a Gaussian with a width that corresponds to the total 
systematic uncertainty. The dominant uncertainties arise from 
the uncertainties on the top quark mass and on the jet energy 
scale (JES). The mass of the top quark is varied by ±2.3 GeV 
and the JES by ±la around their nominal values. 

The statistical uncertainty on the cos 9* templates is taken 
as a systematic uncertainty estimated by fluctuating the tem- 
plates according to their statistical uncertainty, and noting the 
RMS of the resulting distribution when fitting to the data. 

The effect of gluon radiation in the modeling of tt events 
is studied with an alternate MC sample that includes tt events 
generated with an additional hard parton by ALPGEN. These 
events are mixed with the standard tt events according to the 
ratio of the leading order cross sections for these two pro- 
cesses. Effects of the chosen factorization scale Q in the 
generation of the T'F+jets events are evaluated using a sam- 
ple generated with a different choice of Q. The systematic 
uncertainty on the jet flavor composition in the TJ^+jets back- 
ground is derived using alternate MC samples in which the 
fraction of b and c jets are varied by 20% about the nominal 
value l2(ill . The difference found between the input / + value 
and the reconstructed / + value in ensemble tests is taken as 



TABLE II: Systematic uncertainties on /+ for the two channels and 
for their combination. 



Source 


^+jets Dilepton 


Combined 


Jet energy scale 


0.038 


0.039 


0.038 


Top quark mass 


0.019 


0.028 


0.021 


Template statistics 


0.037 


0.024 


0.028 


tt model 


0.006 


0.018 


0.009 


Background model 


0.007 


0.007 


0.005 


Heavy flavor fraction 0.018 




0.015 


Calibration 


0.018 


0.010 


0.016 


Total 


0.063 


0.059 


0.057 



the systematic uncertainty on the calibration of the analysis. 

The systematic uncertainties are conservatively assumed to 
be fully correlated except for those due to template statistics 
and the calibration of the individual analyses, which are com- 
pletely uncorrelated, and the MC model systematic uncertain- 
ties which are partially correlated. Assuming a fixed value of 
0.7 for /o, we find 

/+ = 0.109 ± 0.094(stat) ± 0.063(syst) (2) 
using £+ jets events, and 

/+ = -0.089 ± 0.154(stat) ± 0.059(syst) (3) 
using dilepton events. Combination of these results yields 

/+ = 0.056 ± 0.080(stat) ± 0.057(syst). (4) 

We also calculate a Bayesian confidence interval (using a flat 
prior distribution which is non-zero only in the physically al- 
lowed region of / + = 0.0 — 0.3) which yields 



/+ < 0.23 @ 95% C.L. 



(5) 



As seen in Fig.^a), there is a deficit of €+jets data events 
in the central region of cos 9*. We estimate the significance of 
this effect by performing a likelihood ratio test to evaluate the 
goodness-of-fit for the best-fit model and find that the proba- 
bility of obtaining a worse fit is 1.3%. We also evaluate the 
goodness of fit for the standard model hypothesis and find a 
fit probability of 0.8% (statistical). Thus we conclude that the 
discrepancy is not statistically significant. We have studied 
the subset of our MC ensemble tests in which the mock data 
has a lower fit probability than the collider data does and find 
that our sensitivity to the value of / + in this subset is the same 
as in the entire set of ensembles. 

In summary, we have measured the fraction of right-handed 
W bosons in tt decays in the £+jets and dilepton channels, 
and find /+ = 0.056 ± 0.080(stat) ± 0.057(syst). This is the 
most precise measurement of / + to date and is consistent with 
the standard model prediction of/ + — 3.6 x 10~ 4 |3]. 

We thank the staffs at Fermilab and collaborating insti- 
tutions, and acknowledge support from the DOE and NSF 
(USA); CEA and CNRS/IN2P3 (France); FASI, Rosatom 
and RFBR (Russia); CAPES, CNPq, FAPERJ, FAPESP and 



7 



FUNDUNESP (Brazil); DAE and DST (India); Colciencias 
(Colombia); CONACyT (Mexico); KRF and KOSEF (Korea); 
CONICET and UBACyT (Argentina); FOM (The Nether- 
lands); PPARC (United Kingdom); MSMT (Czech Republic); 
CRC Program, CFI, NSERC and WestGrid Project (Canada); 
BMBF and DFG (Germany); SFI (Ireland); The Swedish Re- 
search Council (Sweden); Research Corporation; Alexander 
von Humboldt Foundation; and the Marie Curie Program. 



[*] On leave from IEP SAS Kosice, Slovakia. 

[f ] Visitor from Helsinki Institute of Physics, Helsinki, Finland. 

[f] Visitor from Lewis University, Romeoville, IL, USA 

[1] G. L. Kane, G. A. Ladinsky, and C.-P. Yuan, Phys. Rev. D 

45, 124 (1992); R. H. Dalitz and G. R. Goldstein, ibid., 1531; 

C. A. Nelson et al., Phys. Rev. D 56, 5 928 (1997). 
[2] TeVatron Electroweak Working Group, hep-ex/0603039 
[3] M. Fischer et al, Phys. Rev. D 63, 031501(R) (2001). 
[4] K. Fujikawa and A. Yamada, Phys. Rev. D 49, 5890 (1994); P. 

Cho and M. Misiak, Phys. Rev. D 49, 5894 (1994); C. Jessop, 

SLAC-PUB-9610. 
[5] CDF Collaboration, T. Affolder et al, Phys. Rev. Lett. 84, 216 

(2000). 

[6] D0 Collaboration, V. M. Abazov et al, Phys. Lett. B 617, 1 
(2005). 



[7] CDF Collaboration, D. Acosta et al, Phys. Rev. D 71, 031101 
(2005). 

[8] CDF Collaboration, A. Abulencia et al, Phys. Rev. D 73, 
111103 (2006). 

[9] D0 Collaboration, V. M. Abazov et al, Phys. Rev. D 72, 
011104(R) (2005). 

[10] D0 Collaboration, V. M. Abazov et al, Nucl. Instrum. Meth- 
ods A 565, 463 (2006). 

[11] Rapidity y and pseudorapidity r\ are defined as functions of the 
polar angle 9 with respect to the proton beam and the parameter 
f3asy(9,/3) = § ln[(l + /3cos0)/(l - (3 cos 9)} and r](6) = 
y(9, 1), where f3 is the ratio of a particle's momentum to its 
energy. 

[12] M. L. Mangano et al, JHEP 07, 001 (2003). 

[13] T. Sjostrand et al, Comp. Phys. Commun. 135, 238 (2001). 

[14] C.A. Nelson et al, Phys. Rev. D 56, 5928 (1997). 

[15] D0 Collaboration, V. M. Abazov et al, Phys. Lett. B 626, 45 

(2005). 

[16] G.C. Blazey et al, hep-ex/00050l5] 

[17] V. Barger et al, Phys. Rev. D 48, 3953 (1993). 

[18] D0 Collaboration, B. Abbott et al, Phys. Rev. D 58, 052001 

(1998). 

[19] D0 Collaboration, V. M. Abazov et al, "Measurement of the 
ti Production Cross Section in pp Collisions Using Dilepton 
Events", in preparation. 

[20] D0 Collaboration, V. M. Abazov et al, Phys. Lett. B 626, 35 
(2005) 



